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Abstract— In this era of modernization, deregulation and 
competition the scenario of the optimization techniques is being 
evolved dramatically towards the functional mimicry of nature. 
There are many nature inspired optimization algorithms such as 
PSO, Ant Colony, Genetic Algorithm, Evolutionary Techniques 
ete. Recently Bacterial foraging Optimization Algorithm has 
attracted a lot of attention as a high performance optimizer. In 
2002, K. M. Passino proposed Bacterial Foraging Optimization 
Algorithm (BFOA)[1]| for distributed optimization and control. 
One of the major driving forces of BFOA is the chemotactic 
movement of a virtual bacterium that models a trial solution of 
the optimization problem. The underlying biology behind the 
foraging strategy of E.coli is emulated in an extraordinary 
manner and used as a simple optimization algorithm. This paper 
presents the BFOA for global optimization. 


Keywords Bacterial Foraging, —_ global 
Chemotaxis, Optimization Algorithm. 
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1. INTRODUCTION 


Optimization is a computational science that studies 
techniques for searching the ‘best’ solutions. It has been 
widely employed in a large variety of fields, including 
transportation, manufacturing, physics, and medicine. To 
tackle complex search problems of the real world, scientists 
have been drawing inspiration from nature and natural 
creatures for years. Optimization is at the heart of many 
natural processes like Darwinian evolution, group behavior of 
social insects, and the foraging strategy of other microbial 
creatures. Natural selection tends to eliminate species with 
poor foraging strategies and favor the propagation of genes of 
species with successful foraging behavior since they are more 
likely to enjoy reproductive success. 

Since a foraging organism or animal takes necessary action 
to maximize the energy intake per unit time spent for foraging, 
considering all the constraints presented by its own 
physiology such as sensing and cognitive | capabilities, 
environment (e.g., density of prey, risks from predators, 
physical characteristics of the search space), the natural 
foraging strategy can lead to optimization and essentially this 
idea can be applied to solve real-world optimization problems. 
Based on this concept, Passino proposed an optimization 
technique known as the bacterial foraging optimization 
algorithm (BFOA) [1]. To date, BFOA has successfully been 
applied to real-world problems such as optimal controller 
design, harmonic estimation, transmission loss reduction, 
active power filter synthesis, and learning of artificial neural 
networks. 
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Il. THE BACTERIA FORAGING OPTIMIZATION ALGORITHM 


In the process of foraging, E. coli bacteria undergo four 
Stages, namely, chemotaxis, swarming, reproduction, and 
elimination and dispersal. In search space, BFOA seek 
optimum value through the chemotaxis of bacteria, and realize 
quorum sensing via assemble function between bacterial, and 
satisfy the evolution rule of the survival of the fittest make use 
of reproduction operation, and use elimination-dispersal 
mechanism to avoiding falling into premature convergence. 


A, Chemotaxis 

This process simulates the movement of an E.coli cell 
through swimming and tumbling via flagella. Biologically, an 
E.coli bacterium can move in two different ways. It can swim 
for a period of time in the same direction, or it may tumble, 
and alternate between these two modes of operation for the 
entire lifetime. Suppose i(j, k, 1) represents i th bacterium at j 
th chemotactic, k reproductive and I" elimination dispersal 
step. C(i ) is the size of the step taken in the random direction 
specified by the tumble (run length unit), Then in 
computational chemotaxis the movement of the bacterium 
may be represented by 


OCF +LAD =O (j,k) +C(i) 


where A indicates a vector in the random direction whose 
elements lie in [-1, 1]. 


B. Swarming 

An interesting group behavior has been observed for 
several motile species of bacteria including E.coli and S. 
typhimurium. When a group of E. coli cells is placed in the 
center of a semisolid agar with a single nutrient chemo- 
effector, they move out from the center in a traveling ring of 
cells by moving up the nutrient gradient created by 
consumption of the nutrient by the group. To achieve this, 
function to model the cell-to-cell signaling via an attractant 
and a repellan. The mathematical representation for E.coli 
swarming [2] can be represented by: 


J(8,P(j.k.)) = ¥ J’ (0,6'(j,k,D) 


i=l 
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where is the cost function value to be added to the actual 
cost function. S is the total number of bacteria and p is the 
number of parameters to be optimized which are present in 


each bacterium. “is the depth of the attractant released 


Ww. i A 

by the cell and = “*' is a measure of the width of the 
h ee 3 

attractant signal. “?*/#"" =“ auract is the height of the repellant 


W,,, - 
effect and = "*?"“”" is a measure of the width of the repellant. 


C. Reproduction 
According to the rules of evolution, individual will 
reproduce themselves in appropriate conditions in a certain 
way. For bacterial, a reproduction step takes place after all 
chemotactic steps. 
N+ 


Jam = YI i kD 
j=l 


Where eee is the health of bacterium i. Sort bacteria and 


chemotactic parameters C(/) in order of ascending cost (higher 
cost means lower health). For keep a constant population size, 


bacteria with the highest J 


bacteria are allowed to split into two bacteria in the same 
place. 


health Values die. The remaining 


D. Elimination-Dispersal 

In the evolutionary process, elimination and dispersal 
events can occur such that bacteria in a region are killed or a 
group is dispersed into a new part of the environment due to 
some influence. They have the effect of possibly destroying 


chemotactic progress, but they also have the effect of assisting 
in chemotaxis, since dispersal may place bacteria near good 


food sources. From the evolutionary point of view, 
elimination and dispersal was used to guarantees diversity of 
individuals and to strengthen the ability of global optimization. 
In BFOA, bacteria are eliminated with a probability of ped.In 
order to keeping the number of bacteria in the population 
constant, if a bacterium is eliminated, simply disperse one to a 
random location on the optimization domain. 


II. THE BFOA ALGORITHM 
[Step 1] Initialization 
i. Number of parameters (p) to be optimized. 

ii. Number of bacteria (S) to be used for searching the 
total region. 

ili, | Swimming length N, after which tumbling of bacteria 
will be undertaken in a chemotactic loop. 

iv. N, the number of iterations to be undertaken in a 
chemotactic loop. (N.>N,). 


Vv. N,- the maximum number of reproduction to be 
undertaken. 
vi. Neg the maximum number of elimination and 


dispersal events to be imposed over bacteria. 
vii. P.g the probability with which the elimination and 
dispersal events to be imposed over bacteria. 


viii The location of each bacterium P(1-p, I-s, 1) which 
is specified be random numbers on [-1, 1]. 
ix. The value of C(i) which is assumed to be constant. 
x, The values of A rrr : Warract . Va anattoe and 


Wyepeltant 


[Step 2] Iterative algorithm for optimization 
This section models the bacteria population chemotaxis, 
swarming, reproduction, elimination and dispersal. 
Elimination-dispersal loop: /=/+/. 
Reproduction loop: k=k+/. 
Chemotaxis loop: j=/+/. 

a) For i=1, 2, ..., S, calculate cost function value for each 
bacterium i as follows. 
= Compute value of cost function J(/, j, k. J). Let Jon( i, 


LAD= Mik D+ TO GK PGA) 


(i.e., add on cell-to-cell attractant effect for swarming 
behaviour). 
© Let SJias= Joti, j, k& Y to save this value since we 
may find a better cost via a run. 
= End for stop loop. 


b) For i=1, 2, ....., S take the tumbling/swimming decision 
¢ Tumble: Generate a random vector A(i)€ KR? with 
each element A,(i) m= 1, 2, ....p, a random 


number on [-1, 1]. 
© Move: let 


6'(j +1,k,) =6',k, + C—O 


Va" (DAG) | 
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Fixed step size in the direction of tumble for bacterium i is 
considered. 


© Compute J(i, 7+1,4,/) and then let Ju( i,j, kD 
= Mii ke D+ S(O kD. PU.KD) 


¢ =Swim: 
[Let m=0; (counter for swim length) 
II. While m<N, (have not climbed down too long) 
e =Let m=m+1 
© fi i+1, kK Y<Sias (if doing better), let jas: Joli, 
j+1, k, )) and let 


Sea 
lntahzanonof 
Vv, 
puss Compute value of 
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Ekmunation and bacterum as J(LP) 
Dispersal Loop where lis the 
counter, L=1 bactenum number 
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Fig. 2(a) Flow chart of Bacterial Foraging Algorithm 
Nutrient concentration (valleys=food, peaks=noxious) 


Fig.2(b) Contour of the test Function 
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A(i) 

VA" (AC) 
and use this 6'(j +1,k,/) to compute the new 
J(i, 7 +1, kD) 

e Else, let m=Ns. This is the end of the while statement. 


c) Go to the next bacterium (i+1) if i # S (i.e. go to b) 
to process the next bacterium. 


OG+LAD =O (Fk, D+ Cli) 


If 7 < N,, go to step 3. In this case, continue chemotaxis 


since the life of the bacteria s not over. 
Reproduction 
a) For the given k and /, and for each i=1, 2, ....S, let 


of = ] 7 sk 
San = amin AS (is j,k, )} be the health of the 
bacterium i (a measure of how many nutrients got 


over its life time and how successful it was at 
avoiding noxious substance). Sort bacteria in order of 


ascending cost J,,..,,, (higher cost means lower 
health). 

b) The S,=S/2 bacteria with the highest J 
die and other S, bacteria with the best value split (and 


the copies that are made are placed at the same 
location as their parent) 


health Values 


If k<N,, go to 2, in this case we have not reached the 


number of specified reproduction steps, so we start the 
next generation in the chemotactic loop. 
Elimination-dispersal : For i=1, 2, ..S, with probability P.g, 
eliminate and disperse each bacterium(this keeps the 
number of bacteria in the population constant) to a 
random location on the optimization domain. 


Flow Chart of the above algorithm is shown in Fig. 2(a) 


IV. SIMULATION AND RESULTS 

The algorithm described above is simulated for the test 
function : 

J = 5*exp(-0.1*((theta(1,1)-15)°2+(theta(2,1)-20)*2))- 
8*exp(-0.08*((theta(1,1)-20)°2+(theta(2,1)-15)*2)) in 
MATLAB 7.9.0. The following results are obtained when 
S=50, p=2, Ne=100, Ns=4, Nre=4, Sr=S/2, Ned=2, Ped=0.25 

Figure 2 shows the contour of the function where valleys 
are food and peaks are noxious. Figure 3 shows the bacteria 
movement towards the global optimization. Figure 4 shows 
the nutrients obtained by bacteria during life. Figure 5 shows 
the health of the bacteria ie. Jhealth . Figure 6 shows the 
location of the bacteria P(1-p,1-s,1). 


V. CONCLUSION 
Bacterial Foraging Algorithm can be effectively used to solve 
optimization problems. The effectiveness of the BFOA is 
verified through simulation results on a test function and # 
provides a scope for more future work on this technique. 
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Fig.4: Nutrients obtain during life 
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